February 2, 2008 16:16 WSPC - Proceedings Trim Size: 9in x 6in WSCarforaBuchert arXivl 



1 



> 



o 



X 



RICCI FLOW DEFORMATION 



00 

o 

O 
(N 

OF COSMOLOGICAL INITIAL DATA SETS 

^ , 

; M. CARFORA* and T. BUCHERTt 

^ Dipartimento di Fisica Nucleare e Teorica, Universita' degli Studi di Pavia, 

^ |. Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 

I . Pavia, Via Bassi 6, 27100, Italy 

i-S^ ' E-mail: mauro.carfora@pv.infn.it 

^ ■ 

T Universite Lyon 1, Centre de Recherche Astrophysique de Lyon, CNRS UMR 5574, 
9 avenue Charles Andre, F-69230 Saint-Genis-Laval, France 
E-mail: buchert@obs.univ-lyonl.fr 



Ricci flow deformation of cosmological initial data sets in general relativity is 
a technique for generating families of initial data sets which potentially would 



allow to interpolate between distinct spacetimes. This idea has been around 



since the appearance of the Ricci flow on the scene, but it has been difiicult 
to turn it into a sound mathematical procedure. In this expository talk we 
illustrate, how Perelman's recent results in Ricci flow theory can considerably 
' improve on such a situation. From a physical point of view this analysis can be 

QQ ' related to the issue of finding a constant— curvature template spacetime for the 

^ inhomogeneous Universe, relevant to the interpretation of observational data 

and, hence, bears relevance to the dark energy and dark matter debates. These 
techniques provide control on curvature fluctuations (intrinsic backreaction 
terms) in their relation to the averaged matter distribution. 



^ Keywords: Ricci flow, Relativistic cosmology, initial-value problem in GR. 



1. INTRODUCTION 

The Ricci flow has been introduced by R. Hamilton^^ with the goal of pro- 
viding an analytic approach to Thurston's geometrization conjecture for 
three-manifolds'^'^''^^ . Inspired by the theory of harmonic maps, he consid- 
ered the geometric evolution equation obtained when one evolves a Rieman- 
nian metric gab, on a three-manifold S, in the direction of its Ricci tensor^'' 
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TZab, i-e. 

(1) 

gabiP = 0)=gab , < /3 < To . 

In recent years, this geometric flow has gained extreme popularity thanks 

to the rcvohitionary breakthroughs of G. Pcrclman^''"^^ , who, taking the 
whole subject by storm, has brought to completion Hamilton's approach 
to Thurston's conjecture. The prominent themes recurring in Hamilton's 
and Perelman's works converge to a proof that the Ricci flow, coupled to 
topological surgery, provides a natural technique for factorizing and uni- 
formizing a three-dimensional Riemannian manifold into locally ho- 

mogeneous geometries. This is a result of vast potential use also in theo- 
retical physics, where the Ricci flow often appears in disguise as a natural 
real-space renormalization group flow. Non-linear u-model theory, describ- 
ing quantum strings propagating in a background spacetime, affords the 
standard case study in such a setting-^'^'^'-^^'^^'^^ . 

Another paradigmatical, perhaps even more direct, application occurs 
in relativistic cosmology^' , (for a series of recent results see also^''' and 
the references cited therein). This will be related to the main topic of this 
talk, and to motivate our interest in it, let us recall that homogeneous and 
isotropic solutions of Einstein's laws of gravity (the Friedman Lemaitre- 
Robertson- Walker (FLRW) spacetimes) do not account for inhomogeneities 
in the Universe. The question whether they do on average is an issue"^^ 
that is the subject of considerable debate especially in the recent literature 
(see^^'^^ and follow-up references; comprehensive lists may be found in^^'^^ 
and^). 

In any case, a member of the family of FLRW cosmologies (the so-called 
concordance model that is characterized by a dominating cosmological con- 
stant in a spatially flat universe model) provides a successful fitting model to 
a large number of observational data, and the generally held view is that the 
spatial sections of FLRW spacetimes indeed describe the physical Universe 
on a sufficiently large averaging scale. This raises an interesting problem in 
mathematical cosmology: devise a way to explicitly construct a constant- 
curvature metric out of a scale-dependent inhomogeneous distribution of 
matter and spatial curvature. It is in such a framework that one makes the 
basic observation that the Ricci flow (1) and its linearization, provide a 
natural technique^'^'^ for deforming, and under suitable conditions smooth- 
ing, the geometrical part of scale dependent cosmological initial data sets. 
Moreover, by taking advantage of some elementary aspects of Perelman's 
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results, this technique also provides a natural and unique way for deform- 
ing, along the Ricci flow, the matter distribution. The expectation is that 
in this way we can define a deformation of cosmological initial data sets 
into a one-parameter family of initial data whose time evolution, along the 
evolutive part of Einstein's equations, describe the Ricci flow deformation 
of a cosmological spacetime. 

2. The kinematical set— up: Initial data set for cosmological 
spacetimes 

To set notation, we emphasize that throughout the paper we shall consider 
a smooth three-dimensional manifold E, which we assume to be closed and 
without boundary. We let C°°(S, M) and C°°(S, (»p r*I](g)«rS) be the space 
of smooth functions and of smooth {p, (j')-tensor fields on S, respectively. 

We shall denote by P?'//(E) the group of smooth diffcomorphisms of E, 
and by 7?,zem(E) the space of all smooth Riemannian metrics over E. The 
tangent space , 7(5],g)^iem(E), to TZiem{T,) at (E, g) can be naturally iden- 
tified with the space of symmetric bilinear forms C°°(E,®^r*E) over E. 
The hypothesis of smoothness has been made for simplicity. Results similar 
to those described below, can be obtained for initial data sets with finite 
Holder or Sobolev differentiability. In such a framework, let us recall that 
a collection of fields g G niem(T.), K e T(^^,g)niem,{j:), g e C°°(S,M+), 
J S C°°(E, TE), defined over the three-manifold E, characterizes a set 
( E , gab , Kab , Q , Ja )j of physical cosmological initial data for Einstein 
equations if and only if the matter fields (g, J) verify the weak energy condi- 
tion (? > 0, the dominant energy condition > gabJ"'J^i and their coupling 
with the geometric fields {g, K) is such as to satisfy the Hamiltonian and 
divergence constraints^: 



Here A is the cosmological constant, k := g"'^Kab, and TZ is the scalar 
curvature of the Riemannian metric gab- If such a set of admissible data 
is propagated according to the evolutive part of Einstein's equations, then 
the symmetric tensor field Kab can be interpreted as the extrinsic curva- 
ture and k as the mean curvature of the embedding it : T, ^ M^^^ of 



''Latin indices run through 1, 2, 3; we adopt the summation convention. The nabla oper- 
ator denotes covciriant derivative with respect to the 3-metric. The units are such that 
c= 1. 



n + k^- K\K\ = IQ-kGq + 2A , 
WbK\ - Vak = SnGJa • 



(2) 
(3) 
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(E, (/ab) in the spacetime (M^^^ ~ S x M, resulting from the evolution 
of (E, Qab, Kab, Q, Ja)i whereas q and Jo are, respectively, identified with the 
mass density and the momentum density of the material self-gravitating 
sources on (E,£fa6)- 

3. The Heuristics of averaging: Deformation of 
cosmological initial data sets 

The averaging procedure described in^'^'i" is based on a smooth deforma- 
tion of the physical initial data {T,, gab, Kab, g, Ja) into a one-parameter 
family of initial data sets 

/3 ^ (E , gabiP) , KabiP) , QiP) , Jam), (4) 

with < /3 < 00 being a parameter characterizing the averaging scale. The 
general idea is to construct the flow (4) in such a way as to represent, as 
f3 increases, a scale-dcpcndcnt averaging of (E, gab, Kab, Q, Jo), and - under 
suitable hypotheses - reducing it to a constant-curvature initial data set 

, Qab > Kab = ^9ab^ ,Q, Ja = 0), (5) 

where is a constant curvature metric on E, A: is the (spatially constant) 
trace of the extrinsic curvature (related to the Hubble parameter), and 'g 
is the averaged matter density. Under the heading of such a general strat- 
egy it is easy to figure out the reasons for an important role played by the 
Ricci flow and its linearization. As we shall recall shortly, they arc natural 
geometrical flows always deflning a non-trivial deformation of the metric 
gab and of the extrinsic curvature Kab- Moreover, when global, they posses 
remarkable smoothing properties. For instance, if, for /3 = 0, the scalar 
curvature TZ of {T,,gab) is > 0, and if there exist positive constants ai, 
a2, as, not depending on (3, such that Ti-abiP) — C(i9ab{f3)T^{P) > 0, and 
Tiab{m''\P) < a^n^-'^'iP), where UabiP) = 7^a6(/3) - yabimiP) is 
the trace-free part of the Ricci tensor, then'^^ the solutions (gab{f3), KabiP)) 
of the (volume-normalized) Ricci flow and its linearization exist for all 
/3 > 0, and the pair {gab{0),Kab{0)), uniformly converges, when /3 — > oo, 
to iVabi ^v9ab) where is a metric with constant positive sectional cur- 
vature, and V is some vector field on E, possibly depending on /?. The flow 
P 1-^ iVab^ ^vdab)' describing a motion by diffeomorphisms over a constant 
curvature manifold, can be thought of as representing the smoothing of the 
geometrical part of an initial data set (E, gab, Kab, Q,Ja)- 
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It is useful to keep in mind what we can expect and what we cannot expect 
out of such a Ricci-flow deformation of cosmological initial data set. Let us 
start by remarking that the family of data (4) will correspond to the initial 
data for physical spacetimes iff the constraints (2) and (3) hold throughout 
the /3 dependent deformation. This is a very strong requirement and, if we 
have a technically consistent way j3 ^ (QiP), gabiP), Kab{P)) of deforming 
the matter distribution and the geometrical data, then the most natural way 
of implementing the constraints is to use them to define scale-dependent 
backreaction fields /? i-^ </'(/?), /3 ''PaiP) describing the non-linear inter- 
action between matter averaging and geometrical averaging, i.e., 



m = q{P) - (167rG)-i [7^(/3) + k\p) - K\{I3)K\{(3) - 2A] , (6) 



To illustrate how this strategy works, let us concentrate, in this talk, on the 
characterization of the scalar field (3 ^ (t>{f3), providing the backreaction 
between matter and geometrical averaging. The covector field f3 i-> tpaiP) 
can, in principle, be controlled by the action of a /3-dependent diffeomor- 
phism. However, its analysis requires a subtle interplay with the kinematics 
of spacetime foliations^ , {i.e., how we deal with the lapse function and with 
the shift vector field in the framework of Perelman's approach), and will be 
discussed elsewhere, (for a pre-Perelman approach to this issue see^'^°) . 

Let us start by observing that the matter averaging flow (3 i— *■ g{f3) must 
comply with the preservation of the physical matter content 



and must be explicitly coupled to the scale of geometrical averaging. In 
other words, if, for some fixed /3 > 0, we consider that part of the mat- 
ter distribution g{P) which is localized in a given region B{x,t{(3)) C Tip 
of size t(/3), then we should be able to tell from which localized distribu- 
tion {Qm,B{x,T)), at /? = 0, the selected matter content {Q{fi), B{x,t{(5))) 
has evolved. A natural answer to these requirements is provided by Perel- 
man's backward localization^^ of probability measures on Ricci evolving 
manifolds. The idea is to probe the Ricci flow with a probability mea- 
sure whose dynamics can localize the regions of the manifold (S, g) of ge- 
ometric interest. This is achieved by considering, along the solution QabiP) 
of (1), a /^-dependent mapping f3 i — > /(/?, ) e C°°(I]/3,M), in terms of 
which one constructs on the measure dzu{0) = (47rr(/?))~^ e~^ d^g(^fj), 



(7) 




(8) 



February 2, 2008 16:16 



WSPC - Proceedings Trim Size: 9in x 6in WSCarforaBuchert'arXivl 



6 

where (3 i — > T{j3) G R"*" is a scale parameter chosen in such a way 
as to normaUze dw{j3) according to the so-cahed Perelman's coupling : 
Jj,^ dw{(3) = (47rr(/3))~^ J^,^ dUgt^ij) = 1. It is easily verified that this 
is preserved in form along the Ricci flow (1), if the mapping / and the scale 
parameter t(/3) arc evolved backward in time /3 € (/3*,0) according to the 
coupled flows defined by 

' ^/ = -A,(^)/ - R{P) + |t(/3)-i, /(/?*) = /o 

(9) 

_ ^t(/3) = -1,t(/3*)=to, 

where ^g(fi) is the Laplacian with respect to the metric gab{P), and /o, tq 
are given (final) data. In this connection, note that the equation for / is a 
backward heat equation, and as such the forward evolution /(/? = 0) ^ / is 
an ill-posed problem. A direct way for circumventing such a difliculty is to 
interpret (9) according to the following two-steps prescription: (i) Evolve 
the metric /3 i — > (S, say up to some /?*, according to the Ricci 
flow, (if the flow is global we may let /?* — » cxd); (ii) On the Ricci evolved 
Riemannian manifold (S,^^^(/3*)) so obtained, select a function /(/?*) and 
the corresponding scale parameter t(/3*), and evolve them, backward in /?, 
according to (9). 



4. Ricci— flow deformation of the initial data {T,,gabiKabi o) 

With these preliminary remarks along the way, let us characterize the var- 
ious steps involved in constructing the flow (4). (For ease of exposition, 
we refer to the standard unnormalized flow; volume normalization can be 
enforced by a reparametrization of the deformation parameter). 

Definition 4.1. [ Geometrical Data Deform,a,tion) Given an initial data set 
{T,,gab,Kab,g,Ja) for a cosmological spacetime (M^^^ ~ E x R, (/(^)), the 
Ricci-flow deformation of its geometrical part (S, gat, Kab) is deflned by the 
flow /3 {gab{l3), Kfij{f]j), < (3 < /3* provided by the (weakly) parabolic 
initial value problem (the Ricci flow, proper) 

^9ab{P) = -2TZab{P), 
< gabiP = 0)=gab , 0<P< P* (10) 



lim/3/'/3« [sup^gs < oo > 
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and by its linearization, which, by suitably fixing the action of the dif- 
feomorphism group 2?i//(S), takes the form of the parabolic initial value 
problem^^ 

(11) 

KabiP = 0)= Kab , 0<P<P* , 

where Al denotes the Lichnerowicz-DeRham Laplacian A^Kab = 
V^ViKab — RasK^ — RbsKl + 2RasbtK^*' acting on symmetric bilinear 
forms^^ . 

Definition 4.2. {Localization of the Deformed Data) The geometrical de- 
formation P {gabiP)^ Khj{P)) is controlled by the backward localizing 
flow rj {E{x,y,ri),E°;\,,{x,y,ri)), rj = 0* — l3, defined by the backward 
heat kernels for the conjugate^ ^'^^ parabolic initial value problem associ- 
ated with the g((?7)-dependent Laplace-Beltrami and Lichnerowicz operators 
77 1-^ (A, Ai). The operator A^, when acting on (bi)scalars, reduces to A. 
Thus, we can characterize both these kernels in a compact form as solutions 
of 

r ^ (EfAx,y,r^)\^ (X) fE-,i,{x,y;v)\ _ ^ (E<l\,{x,y;r,)\ 
V E{x,y;n) ) ^ \ E{x,y;v) J \ E{x,y;v) J ' 

< 

fEf,^,{x,y;r,)\ _ f6f,^,ix,y;rj)\ 

(12) 

where {y,x;r]) G (S x j:\Diag{T, x S)) x [0,/3*], A^""^ denotes 
the Lichnerowicz-DeRham Laplacian with respect to the variable x, 
the heat kernels E°;\,{y, x;ri) and E{x,y;ri) are smooth sections of 
((gj^TS) Kl ((g)2T*S) and EKE, respectively, and finally, x; 77) 

is the Dirac p-tensorial measure, {p = 0,1,...) on {T,,g{r])). The 
Dirac initial condition is understood in the distributional sense, i.e., 
/e -^i'fe' ^! ^) ^' '^'(y) ^l^^g{ri) ~* w°'''{x) as r? \ 0+, for any smooth 
symmetric bilinear form with compact support w;* <E C^(E, (g^TE), and 
where the limit is meant in the uniform norm on C^(Tt, ®^TE). 

Note that heat kernels for generalized Laplacians, such as A^, (smoothly) 
depending on a one-parameter family of metrics e 1-^ dabi^), e > 0, are 
briefly dealt with in^ . The delicate setting where the parameter dependence 
is, as in our case, identified with the parabolic time driving the diffusion of 
the kernel, is discussed in^^'^^ , (see Appendix A, §7 for a characterization 
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of the paramctrix of the heat kernel in such a case), and in^^ . Strictly 
speaking, in all these works, the analysis is confined to the scalar Laplacian, 
possibly with a potential term, but the theory readily extends to generalized 
Laplacians, always under the assumption that the metric QabiP) is smooth 
as y (3* ■ Finally, the kernels for A and can both be normalized, along 
the Ricci flow, over the romid (collapsing) 3 sphere (S'^, .9can(/3))• 
Let us now consider the matter content localized by taking the dwirf)- 
expectation of q{iii), 

J^Q{rj)dvuiri)=M{dvu{rj)) . (13) 

According to (8), we require that such a local mass is preserved along the 
ry-evolution of the measure dw{r]), i.e., Jj, Q{r]) dw{r]) = 0. This request 
motivates the following 

Definition 4.3. {Deformation of Matter Data) The given matter distri- 
bution g{(3 = 0) is deformed according to the heat-flow (3 i— > q[(3) given 
by 

^K/3) = Ag(^) g{f3), /3e[0,To) 

(14) 



We are now in the position to characterize the Ricci flow deformation of 
the cosmological data (S, gab, Kab, Q, Ja)- According to the results described 
in^^ we have 

Proposition 4.1. Let /3 i— > {gab{f3),Kik{/3),Q{(3)) be the Ricci flow defor- 
mation, on S,, X [0, P*], of the data (S, gab, Kab, q) as defined above. Assume 
that the underlying Ricci flow [3 > (S, gabiP)) is of bounded geometry, and 
let E{y, x; rj) and E^il],, {y, x; rf) be the (backward) heat kernels Ricci-flow con- 
jugated to the Laplace-Beltrami and to the Lichnerowicz-DeRham operator, 
respectively. Then, for all <ri < p*, 

9i'k' {y,V = Q) = Efk' {y, x; rf) [gab{x, r?) - 2// Hab{x, r?)] rfMg(x,r,) , (15) 

KMy,v = o)j = y,^-'(^'-;^) [Kabix,,)) '^^^^-^ ' 

and g{y,r] = 0) ^ j E{y,x;vi) q{x,71i) djXg^^^^) . (17) 
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Moreover, as r] \0~^, we have the uniform asymptotic expansion 



^ [ exp (_^) (l^^f^l) a,,,.,,, 

(47r77)5 Je V 4?7 7 \Kab{x,r])J 

S L J ^t'^-' U.b(., j '^^^^-^ 



+0(77^-^) , (18) 

w/iere t°;\, {y , x; rj) G TS^ KIT*E^ zs the parallel transport operator associ- 
ated with (T,,g{ri)), doly^x) is the distance function in {T,,g{rj = 0)), and 
^h]f^,{y,x;ri) are smooth sections G C°°{Y. x S', ^^TS H (g)2r*S), (de- 
pending on the geometry of {T,, g{r]))), characterizing the asymptotics of the 
heat kernel E^\,{y,x;r]). With an obvious adaptation, such an asymptotic 
behavior also extends to git^i {y,r] = 0) and g{y,ri = 0). 

With these results it is rather straightforward to provide useful characteri- 
zations of the (intrinsic) backrcaction field (t>{(3). For ilhistrativc purposes, 
let us assume that (3 t-^ Kab{(3) = and A = , then from (6) and the 
above proposition wc get that 

(j){y,r} = Q)= / E{y,x;ri)Q{x,rf}dHg(^^^) (19) 



-(167rG)-i 



^ 5*'"' {y, = 0) Ef^,, {y, x; ry) 7^„b(x, t?) dfx 



If we further assume that the Hamiltonian constraint holds at the 
fixed observative scale r] and that TZabix, rj) w WlZirf) gab{x,r] = 0) + 
5Tiab{x,rj)], [i.e., curvature fluctuates around the constant curvature back- 
ground TZi'm' iv, rj = Q) = 2C gv^' iv, '7 = 0)), then 

</)(y,r? = 0)«(167rG)-iy g"^'"' {y,r] = 0) Efk,{y,x;v) Snab{x,v) dfig^^^r,) , 

(20) 

from which it follows that the intrinsic backrcaction field <^(/3) is generated 
by curvature fluctuations around the given background, as expected. 
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